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Hydrodynamical analysis of experimental data of ultrarelativistic heavy ion collisions seems to 
indicate that the hot QCD matter created in the collisions thermalizes very quickly. Theoretically, 
we have no idea why this should be true. In this proceeding, I will describe how the thermal- 
ization takes place in the most theoretically clean limit - that of large nuclei at asymptotically 
high energy per nucleon, where the system is described by weak-coupling QCD. In this limit, 
plasma instabilities dominate the dynamics from immediately after the collision until well after 
the plasma becomes nearly in equilibrium at time t ^ ^^^Q^^ ■ 
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1. Introduction 

The hydrodynamical description of heavy ion collisions at RHIC and at LHC has led to a 
tremendous phenomenological success. However, the hydrodynamical treatment can be justified 
only if the matter created in the collision is near local thermal equilibrium (or more precisely close 
to local isotropy). At the initial stages of the collision, this condition is clearly violated and it is 
an open theoretical question, how quickly — and how — the matter approaches equilibrium. It is 
crucial to understand this prethermal evolution as our ignorance of it constitutes one of the largest 
systematic uncertainties in analysis of the heavy-ion collisions. 

In [|l|], we have addressed this problem in the most theoretically clean limit — that of large 
nuclei at asymptotically high energy per nucleon, where the system is described by weak-coupling 
QCD. When the typical energy scale right after the collisions is large ^ Aqcd, the renormalized 
strong coupling constant becomes small a (2s) ^ 1. In this hmit, the initial condition for the col- 
lision is understood in the color glass condensate (CGC) framework [^. In |[l]], we have identified 
the most important physical processes, in the a ^ 1 limit, that drive the evolution from the initial 
strong and anisotropic CGC fields to the thermal state in the longitudinally expanding geometry of 
a heavy-ion collision. Our solution resembles the original "Bottom-Up" themialization [^, with 
the difference that it takes into account the physics of plasma instabilities 

The evolution proceeds in three stages. The first stage (1^2sT^a^7)is characterized by 
strong fields, or equivalently high occupancies (/ » 1). It is a competition between the longitudinal 
expansion that drives the system towards larger anisotropics ^ (jPzl) = ol''{\pi_\)) and 

weaker fields, and momentum broadening due to interactions that works towards isotropizing the 
fields. The result of this competition is that the anisotropy increases and the occupancies decrease 
as a function of time. 
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During the second stage <^QsT <^ , the system is highly anisotropic but the typical 
modes are now under-occupied, / ^ 1. The cross-over from high to low occupancy changes quali- 
tatively the system's behavior: inelastic scattering begins to increase particle number. In particular, 
soft inelastic emissions create a bath soft gluons that eventually becomes nearly thermal. During 
this stage, the soft protothemial bath does not dominate anything; the primordial hard particles 
carry the most energy, inflict most screening, are more numerous, and cause the most scattering. 

Eventually the soft protothermal bath will, however, start to dominate the physics, and in 
particular the momentum broadening experienced by the hard particles. This is the third stage 
Q,-i2/5 ^ Q^j. ^ o;^5/2 Qyj^^g this Stage the hard anisotropic particles still carry most of the 
energy density of the system (and hence dominate T^y, relevant for hydrodymamics), but they 
effectively decouple from each other as the dominant interaction is with the bath of soft particles; 
the hard particles can be seen as few but highly energetic "jets" propagating in a nearly thermal 
medium. The thermalization then proceeds through quenching of these jets: interaction with the 
protothermal medium leads to hard coUinear splitting and therefore to radiative energy loss. Once 
the hard modes have had time to lose all their energy to the medium (by the time QgT ~ a^^/^), all 
that is left is the nearly thermal bath of soft modes, and the system has essentially thermalized. 

Throughout the evolution, the dominant interaction between the hard particles, and between 
the hard particles and the protothermal bath, takes place via plasma (or Chromo-Weibel) unstable 
modes: in an anisotropic system, a set of long-wavelength chromo-magnetic fields is perturbatively 
unstable and the unstable modes undergo an exponential growth until they become non-perturbative 
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(with /(^inst) ~ l/a) and saturate. When hard particles propagate thi^ough a stochastic background 
of these saturated 5-fields, they exchange momentum with the B-fields due to Lorentz force. As 
the B-fields are at long scales incoherent, this momentum transfer is diffusive and is described by a 
(time dependent) momentum diffusion coefficient ^i^st ^ dp^ /dt. This q turns out to be larger than 
that arising from ordinary elastic scattering at all stages of thermalization. 

During the first and second stages, the most important instabilities are those due to the hard par- 
ticles. The soft protothermal bath is also sightly anisotropic (because of expansion and anisotropic 
"rain" of soft particles arising from splitting of the anisotropic hard particles), and during the third 
stage it is the magnetic fields that became unstable because of the protothermal bath that dominate 
the momentum transfer. 

Here, we discuss these three stages in slightly more detail, but with much less detail than [jl|], 
concentrating only on the dominant scales. 

2. Initial condition in heavy ion collisions: Color-Glass-Condensate 

In the weak coupling limit, the very eaiiy dynamics are well understood in the color glass 
condensate (CGC) framework [||]. It indicates that at time 2s 'i' ~ 1> the system consists of intense, 
nearly boost invariant gluon fields, with a coherence length /coh 27 ' in the ;c3;-plane transverse 
to the beam axis, and a much longer coherence length in the z-direction along the beam axis, and 
that the energy density of the system is e{r ~ 27') ~ 2i/o; . 

At times Qs'C> 1, the fields have lost phase coherence, and can be described in temis of particle 
degrees of freedom. The corresponding distribution function of gluons can be parametrized by^ 

f{p)^a-'d{Q,-\p\)d{a'Q,-p,), d>0 (2.1) 

where c describes the typical occupation number, and d parametrizes the anisotropy. In terms of 
these descriptors, the energy density condition implies that on the ccf-plane, at times 2s 1' ^ 1, the 
system lies onad = c — I line, as shown in Fig. |l|. We take this as our initial condition. 

Where precisely on the line the initial condition sits is a matter of taste. Both, our description 
and that within CGC framework [||], displays that the system evolves to ^(1) anisotropy in a time 
T ~ 27' X (logs of a). Therefore, at times zQs » 1, but less than any negative power of a, the 
system has c = 1 and d = 0. 

3. Evolution to small occupancies 1 ^ Q^T ^ a~'^/^ 

As long as c = 1 and d = 0, the dynamics of the system are fully non-perturbative, and 
amendable only to non-perturbative (yet classical) simulations [^, |6|]. However, the system is only 

'The -functions should not be interpreted as sharp step functions but rather as smooth momentum cutoffs so that 
the actual form of the distribution is such that high momentum particles with p > do not dominate anything, i.e. 
f(p > Qs) < [Qs/ p)^. This ansatz for the distribution is accurate enough for parametric estimates (of powers of a), 
for full numerical description one needs to consider the details of the distribution function. Also, it assumes that 
the physics of the system is dominated by a single scale only, which need not to be the case. However, multi-scale 
systems can be constructed by superimposing several distribution functions of this kind, as we will do in Section ^ The 
functional (power-law) form of / is described in detail in here we limit ourselves to this simpler description for 
clarity. 
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Figure 1: (Left) Graphical representation of the descriptors of an anisotropic single-scale system. (Right) 
The CGC initial condition lies of a c/ = c — 1 line on cc/-plane. 

marginally non-perturbative; deviation from this point by any positive power of a (c < 1 or J > 0) 
renders parts of the system perturbative. And indeed, the interactions and the expansion both work 
toward driving the system away from this point. 

Once the system is away from the non-perturbative point, several scale separations emerge. In 
particular, the kinetic mean free path of particles with p ~ becomes much larger than their ther- 
mal wavelength, Ifree ~ o;^^^''+^'^27' » 27^, so that their evolution and mutual interactions can 
be described in an effective kinetic theory [^. The screening scale (ni^ J d^pf{p) //? ~ cc'^^^ Q^) 
becomes parametrically longer than the wavelength of the typical hard particles. The interaction 
between the modes at scales m and 2s is still non-perturbative, but between modes at scale m 
perturbative. The non-perturbative interaction between m and scales can be resummed by de- 
scribing the modes in the screening scale by coUisionless Vlasov equations, or equivalently in the 
hard-loop (HL) effective theory. The non-equilibrium system also has an infrared scale, analogous 
to the magnetic scale (aT) of the thermal ensemble, characteiized by non-perturbative interactions 
between modes at the infrared scale. Fortunately, this scale does not dominate anything.^ 

3.1 Longitudinal expansion 

The effect of the spatial expansion translates into redshift of the p^ components of the momenta 
dp,/dt ~ —Pz/t and reduces particle number n'{t) '-^ —n{t)/t. If the system is highly anisotropic 
d >0, the energy E ~ ^ p^ + p^j_ of the hard paiticles is dominated by pj^, and the red-shifting does 
not appreciably affect the particle energies. Then also the energy density scales as e'{t) ~ —e{t)/t. 
Therefore, at later times energy conservation forces the system to be constrained on other lines of 
fixed (smaller) energy density on the cJ-plane, J = c — 1 + a for Q^T ~ a^" , as displayed in the 
Fig. m Where exactly on these lines the system takes its place depends on the details of dynamics. 

3.2 Momentum broadening and anisotropic screening 

An 1^(1) part of momentum transfer experienced by a hard particle propagating through an 
equilibrium plasma comes from momentum exchanges of the order of the screening scale. In an 
anisotropic plasma there are magnetic modes at the screening scale which become large, and hence 
the soft momentum transfers in an anisotropic system become enhanced and even more important. 

In an isotropic medium with f{p) ~ /iso(|p|), long wavelength chromo-electric fields are sta- 
bilized by the physics of screening: the introduction of a background electric E field deflects 



^In [g] it has been argued that the infrared scale might dominate particle number and contain a Bose-Einstein 
condensate. For further discussion in this topic see [pp. 
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Figure 2: (Left) Solutions to the evolution in ctZ-plane. In absence of interactions the system follows 
"Free streaming" line. Including perturbative scattering but neglecting instabilities leads to the attractor 
labeled "elastic" from Inclusion of the effect of instabilities leads to more isotropic attractor labeled 
"Instabilities". (Right) Cartoon of particle propagating through plasma unstable magnetic fields. Each patch 
of same-sign B-field gives rise to a momentum transfer A/j^ck ^ gBl^dh, leading to ^ (A/5kick)^/'coh- 



trajectories of hai^d particles in a way that induces a current J; the electric field is reduced by 
the current, eventually canceling it. Now, the energy originally deposited in the electric field is 
stored in currents, which start to create an electric field that in its turn will quench the currents. 
Thus, the energy oscillates between a charge separating deformation of the distribution of hard 
particles and the electric field. The frequency of the plasma oscillation is related to the screen- 
ing scale Wpjly;^^ ~ = a / d^pf /p, and the dispersion relation of modes with p ^ m becomes 

^(p) ~ \J + ^pi- isotropic medium, static magnetic fields are not screened: the deforma- 
tion to the hard particle distribution due to a static B-field is simply an overall rotation around the 
axis set by a magnetic field. The deformed distribution is identical to the original one, and therefore 
no net currents are created, and correspondingly the "magnetic plasma frequency" Wmagl/iso ^' 
magnetic fields are neither stabilized nor destabilized (conesponding to w^ag < 0) by screening. 

In an anisotropic plasma with an angle dependent particle distribution /aniso(^) the rotation 
does not leave the hard distribution unchanged and leads to non-zero currents. Therefore the B- 
fields may be stabilized or destabilized. For an anisotropic distribution, we may nevertheless ask 
what is the effect averaged over the direction B = B/\B\ (and polarizations) of the magnetic field, 
/aniso(p))B- "^^^^ is equivalent to angle averaging over the particles' momentum distri- 
butions, and hence corresponds to some isotropic system with /ave = {f{p))p, which is neither 
stabilized nor destabilized ft)^ag(^)l/avc = 0. Therefore, even in anisotropic systems, the medium's 
impact on magnetic fields, averaged over all directions, is neutral.^ However, if there are any di- 
rections that ai^e stabilized (OD^ag > 0), there must be other directions that have the opposite effect 
and are destabilized (ft)^a<j < 0). These are the plasma unstable modes. They are always present in 
anisotropic systems, and they grow exponentially B{t) ~ B(0)e^ with a growth rate yr-^ m. 

Which modes become unstable depends on the details of the anisotropic distribution. The set 
of unstable modes can be found, e.g., by finding the range of momenta for which the retarded HL 
propagator has poles in the upper half complex-plane. This question has been addressed in [|^, pi] ]. 



and for the distribution of Eq. ( 2. 1 ), the unstable modes have k'-^^^ < m and k™^^ ~ a m 



^In terms of HL effective theory, this is to say that the angle average of the HL gluon polarization tensor is propor- 
tional to the isotropic, thermal polarization tensor. 
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Nothing grows exponentially forever, and the growth of the magnetic fields stops when some 
new physics kicks in. The previous discussion relied on linearizing in the size of the magnetic field; 
once dynamics become non-linear, the non-perturbative interactions stop the exponential growth 
and the instabilities saturate. The non-linear physics enters through the covariant derivative = 
dfi + igA^, and therefore the non-linear physics kicks in when the interaction term competes with 
the derivative in D^. In a non-abelian theory, the magnetic modes can interact among themselves, 
so the relevant momentum scale is that of the instabilities ~ k^^^^, corresponding to A ~ k^l^^/g, or 
b2 ^ (A X A)2 ~ {k'"''Af ~ k'l'^ klJa, or equivalently / ~ 1 /a.^ That the non-linear interaction 



indeed saturates the instability has been observed in numerical HL simulations ||12|]. 

How important for dynamics the unstable magnetic fields are depends on how large they grow 
and how strongly correlated they are. When the hard particles move through the large magnetic 
fields, they experience a time varying Lorentz force (see Fig. §). The magnetic field remains 
coherently in the same direction for the characteristic coherence length of the magnetic field, 
^coh ^"l**' for a particle moving in a ± direction. The momentum accumulated in one co- 
herence length is Apkick ~ gSIcoh, and propagating through many such uncorrelated patches of 
coherent magnetic field leads to diffusive momentum transfer described by a momentum diffu- 
sion coefficient ^in^t ~ (Apkick)V/coh ~ aB^/k^l'' ~ ^Lt^±" ~ oc'^''m^ ~ a-^+^-^^-'^Ql As ex- 
pected, this plasma instabihty induced q is larger than the one due to perturbative elastic scattering 
^ei ~ a^^'-'^'^Ql during the equilibration process, and therefore we do not need to discuss elastic 
scattering in the following. 

3.3 Competition between momentum broadening and longitudinal expansion 

We are now ready to solve the trajectory of the system in the cJ-plane during the first and 
second stages. The longitudinal expansion makes the distribution more anisotropic with a rate 
dp- /At = —p,/T corresponding to dt//da| expansion = +1, while the momentum broadening isotropizes 
the distribution with a rate dp^/Ax ~ qinst/ Pz corresponding to dJ/da|inst ~ —a^i'^+ii^^^y" ~ 
+a^'*'^'+2", with the energy conservation condition I — c = a — d. These two opposite effects can 
compete only if also dd/da\inst ^{oc^^), which happens when datt. {a) = |, and Catt. = 1 — \a. This 
is solution plotted on the cJ-plane in Fig. ^.^ 

The solution is attractive in the sense that if the system is in a state below (above) the attractor 
{d < |), the system becomes more anisotropic (isotropic) and reaches the attractor along a line of 
constant energy density in a time scale comparable to the age of the system. 

In Fig. ^ we also plot the con^esponding attractor neglecting the effect of plasma instabilities, 
so that the momentum broadening is due ordinary elastic scattering. As the interaction is weaker 
the distribution is, at a given time, less isotropic than with instabilities. The attractor has a comer 
at / 1 where the system turns from over- to under-occupied and q^i °^ f{^+ f) changes behavior. 
For a plasma instability driven system, there is no such kink in the attractor; the plasma instabilities 
depend on the hard particle distribution only through the scale m, which is linear in /, and the 
attractor is valid also at Q^x » a^7. If no new physics would kick in, the system would not 
thermalize and would only get more and more dilute as a function of time. 



'^In |jl|], a more precise gauge invariant criterion is presented in ternis of Wilson loops. For here, this simpler criterion 
is however good enough. 

^We also find another, weakly anisotropic attractor For details, and why we do not think it is realized, see [|l|]. 
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Figure 3: (Left) Cartoon of the infrared cascade of a particle with momentum k. Once the particle has had 
time to split in two daughters with comparable momenta, the daughters quickly cascade to the infrared until 
they reach the soft thermal bath. (Right) Cartoon of the momentum scales in late stages of the evolution. 
The important scales are (from hard to soft): highly anisotropic distribution of hard particles at scale 
(pink), scale k^^^xii below which particles have had time to cascade to infrared, soft nearly thermal bath with 
anisotropy e (black and pink), and the anisotropic screening scale of the T-bath e^^^m with strong magnetic 
fields (dark red). The system becomes effectively thermalized once ^spiit reaches gs- 

__12 5 

4. Creation of a soft nearly thermal bath, "Bottom-up", a s <^ Q^z <^ a^i 

During the first stage, inelastic scattering in an over-occupied system works towards joining 
the hard particles and cascade energy to the ultraviolet. Once the system becomes under-occupied, 
this behavior changes qualitatively and the inelastic scattering starts to split hard particles and 
cascade their energy to the infrared. This leads to the creation of a two-scale system; in addition to 
the original population of hard particles, the debris of splitting form a new, soft^, population. The 
constituents of the soft bath have less inertia, and hence they may isotropize and thermaUze faster 
(by the time Qsl" ~ ot^^i , see [|I|] for details). 

How much energy is deposited in the soft bath and what its temperature is depend on how 
effectively the hard particles cascade to the infrared. The rate at which a hard particle undergoing 
transverse momentum diffusion emits daughters of momentum pdaught. is given (in the LPM regime) 
'^y ^spiit ~ 01 Y^^/ Pdaught.- That is, by the time T all particles with momentum k < ^spUt ~ a^q'C^ have 
had enough time to emit a daughter whose energy is compai^able to the emitter's, corresponding 
to splitting the original particle democratically into two daughters with half the original energy 
(see Fig. The daughters have a higher splitting rate than the mother particle, and they undergo 
successive democratic splittings in a time that is shorter than the age of the system, cascading their 
energy to the soft bath. As by the time t each hard particle have had time to emit particle with 
energy A:spiit, the energy density and the temperature of the soft bath is e ~ T^^^^ ~ ^spiit / d^pf{p)^ 

Even if the soft sector does not dominate the energy density (/ d^pf{p)p), it can dominate 
other characteristics of the medium, such as screening (/ d^pf{p) / p) and q. When q is dominated 
by interactions with the soft bath, the hard particles effectively decouple from each other and see 
only the soft medium. In this case the physical picture is that of a system consisting of a nearly ther- 
mal bath through which a distribution of few but highly energetic "jets" with energy propagate. 
The interaction with the medium quenches the jets via radiative energy loss. When the jets have 

''Soft compared to but still hard compared to the screening scale 

7 — — 

For Qst ^<X 5 , ^ is a strongly angle dependent function giving rise to subtleties not relevant for discussion here. 
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had time to lose all their energy to the medium, that is ^spiit ~ 2s or equivalently Tgq \/ Qs/q, 

the system has effectively thermalized. 

How fast the system then thermalizes is then controlled by how strongly the soft nearly thermal 
bath broadens the momentum of the jets. In [Q], it was assumed that the primary mechanism is 
through perturbative elastic scattering, so that ~ oi^T^^^^. At the time scale the hard jets deposit 
their energy to the thermal bath, the energy density of the thermal bath is comparable to that of the 
hard jets, T^^^^ ~ OL^^ Qt / {Q&'^) , leading to an estimate for the thermalization time Q^'Cq ~ a^^^l^ . 

However, also in this case the plasma instabilities dominate the momentum transfer. The 
soft nearly thermal bath is also anisotropic, partly due to the "stretching" by expansion and partly 
because the cascade of the hard particles heating the bath is anisotropic. Both of these effects lead 
to a parametrically weak anisotropy of the soft bath, of order e = {\pi_\) / {\Pz\) ~ Iri 
[|l|], we found that weakly anisotropic systems such as this generates unstable modes in a range 
where ~ aT^ is the screening scale of the soft distribution.^ These 
unstable modes then subsequently cause momentum broadening ~ k^^^^ ~ £^/^a^/^r^, which is 
larger than the elastic contribution arising from the soft sector as long as £ 2> a^/^, and dominates 
over q arising from the hard distribution for Q^x ^ a^^^l^ . See Fig|3|for cartoon of the scales. 

Taking everything together and solving self-consistently gives 

^sput ~ «5(2,(e,T)2, T ~ aQ,{Q,zfl\ qe ~ a'Ql e ~ a-\Q,r)-'/\ 

so that the hard particles have had time to cascade to the thermal bath, i.e. ^spiit ~ 2s» by the 
time QsT ~ a^^/^. At this point the system (especially its T^v) is approximately isotropic and is 
amendable to a hydrodynamical description. 

References 

[1] A. Kurkela and G. D. Moore, JHEP 1112 (2011) 044 [arXiv: 1 107.5050 [hep-ph]]; A. Kurkela and 
G. D. Moore, JHEP 1111 (2011) 120 [arXiv: 1108.4684 [hep-ph]]. 

[2] E. lancu and R. Venugopalan, In *Hwa, R.C. (ed.) et al.: Quark gluon plasma* 249-3363 
[hep-ph/0303204]. 

[3] P. Romatschke and R. Venugopalan, Phys. Rev. Lett. 96 (2006) 062302 [hep-ph/05 10121]. 

[4] R. Baier, A. H. Mueller, D. Schiff and D. T. Son, Phys. Lett. B 502 (2001) 51 [hep-ph/0009237]. 

[5] S. Mrowczynski, Phys. Lett. B 214 (1988) 587 [Erramm-ibid. B 656 (2007) 273]; 

[6] J. Berges and S. Schlichting, Phys. Rev. D 87 (2013) 014026 [arXiv: 1209.0817 [hep-ph]]. 

[7] P B. Arnold, G. D. Moore and L. G. Yaffe, JHEP 0301 (2003) 030 [hep-ph/0209353]. 

[8] J. -P Blaizot, F. Gehs, J. -R Liao, L. McLerran and R. Venugopalan, Nucl. Phys. A 873 (2012) 68 
[arXiv: 1107.5296 [hep-ph]]. 

[9] A. Kurkela and G. D. Moore, Phys. Rev. D 86 (2012) 056008 [arXiv: 1207. 1663 [hep-ph]]. 

[10] P Romatschke and M. Strickland, Phys. Rev. D 68 (2003) 036004 [hep-ph/0304092]. 

[11] P B. Arnold and G. D. Moore, Phys. Rev. D 76 (2007) 045009 [arXiv:0706.0490 [hep-ph]]. 

[12] P B. Arnold, G. D. Moore and L. G. Yaffe, Phys. Rev. D 72 (2005) 054003 [hep-ph/0505212]. 

^In fact the fc'"*' is given by the screening scale of the anisotropic part of the distribution irq^^^^ ^ saT^ as the 
isotropic part does not stabilize nor destabilize the modes. 



8 



